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Abstract 

An alternative equilibrium stochastic dynamics for a Brownian particle in inhomogeneous space 
is derived. Such a dynamics can model the motion of a complex molecule in its conformation space 
when in equilibrium with a uniform heat bath. The derivation is done by a simple generalization of 
the formulation due to Zwanzig for a Brownian particle in homogeneous heat bath. We show that 
if the system couples to different number of bath degrees of freedom at different conformations then 
the alternative model is derived. We discuss results of an experiment by Faucheux and Libchaber 
which probably has indicated possible limitation of the Boltzmann distribution as equilibrium 
distribution of a Brownian particle in inhomogeneous space and propose experimental verification 
of the present theory using similar methods. 
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Complex molecules like proteins, colloids etc in contact with a homogeneous heat bath 
form a unique class of statistical systems. The most important characteristic feature of 
such systems is that such systems can be modelled by a Brownian particle (BP) with co¬ 
ordinate dependent damping [l-3|. In different conformations (structural configuration of 
a complex molecule e.g. a protein), the accessible surface area of such a molecule to bath 
degrees of freedom can be different. Coupling of the molecule to the uniform heat bath 
can then be conformation dependent 4^6]. The motion of the complex molecule in its con¬ 
formation space under equilibrium fluctuations is then similar to the motion of a BP in a 
multidimensional space with coordinate dependent damping. The source of inhomogeneity 
being damping, even in the absence of any homogeneity breaking potential (conservative 
force), the conformation space remains inhomogeneous. This is a situation for equilibrium 
statistics where there exists an inhomogeneous space in a homogeneous heat bath and no 
extra energetics is involved in maintaining the inhomogeneity of space. This is different from 
an inhomogeneous heat bath where to maintain the inhomogeneity of the heat bath some 
external agent (presence of at least a third system) must be involved and one cannot talk of 
equilibrium without taking into account the presence of that external agent and the entropy 
produced by that involvement. 

Relevant questions one can ask, given such a system, are - a. If the equilibrium probability 
distribution is the Boltzmann distribution when the damping is constant over space, what 
is the general probability distribution when damping is a function of space? This is an 
important question because the damping term cannot be accommodated in a Hamiltonian, 
b. When damping is inhomogeneous, what is the diffusivity? This is a relevant point 


iffusivity is a non-local quantity 
7]. c. How should one generalize 


because, damping being local can be inhomogeneous, but, d 
and should, therefore, be some average over space and time 
the expression of stochastic force to decouple the local damping from the non-local diffusivity 
in equilibrium? As the present paper will elaborate, in the answer of this third question lie 
the answers of the previous two questions. 

In the conventional approach, coordinate dependent damping makes the equilibrium 
stochastic problem involve multiplicative noise. In so far existing literature on such the¬ 
ories the condition for equilibrium is considered to be the appearance of the Boltzmann 
distribution. This is questionable because the distribution does not involve damping which 
also breaks the homogeneity of space. Moreover, these standard procedures give rise to a 
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host of other issues involving the dilemma of Ito vs Stratonovich conventions BO etc. 

Let us first introduce the alternative approach here that has been introduced and elabo¬ 
rated in ref.jl2, 13]. Consider the Brownian dynamics over a homogeneous space character¬ 
ized by a damping constant L 


dx 

dt 


F{x ) 


fm- 


In the above expression F(x) = is the conservative force and g = \j2YkBT is 

the stochastic noise strength where T is the temperature and fog is the Boltzmann con¬ 
stant. The Gaussian white noise 77 (f) of unit strength has correlations < r)i(t) >= 0 and 
< r}i(ti)r]j(tj ) >= SijS(ti — tj ) as usual. The equilibrium probability distribution for the po¬ 
sition of the BP given as P{x) = IVe - ' ( a )/ fc s T where N is a normalization constant. Taking 
an average on the above equation one gets < 777 >= 0 where the noise term vanishes on 
average and the conservative force induced current identically vanishes on average using the 
Boltzmann distribution. 

Now consider a generalization of the above model where T(x) and g(x) are space depen¬ 
dent functions as 

dx F{x) g(x) 

~dt ~ T(x) + f(x) 

Choosing g{x) = \/2Y{x)k^T and 77 (f) being the same zero average white noise one conven¬ 
tionally obtains exactly the same Boltzmann distribution for the so called equilibrium of the 
particle. Take an average of the above equation again and < F(x)/Y(x) > 7 ^ 0 is the case 
now. The average over the noise term is then supposed to cancel out this non-zero current 
contribution in equilibrium. 

The noise term is locally Gaussian with a locally fixed width and zero average. A system 
to equilibrate over an inhomogeneous space like the given one the space must be of finite 
extent even in the absence of the conservative force. Given a long enough trajectory which 
revisits a particular coordinate many times, the noise term in principle can locally everywhere 
sample enough realizations to become increasingly small on average everywhere as the length 
of the trajectory increases. This simple fact means that given a locally Gaussian noise 
everywhere the noise term can tend towards zero on average everywhere over increasingly 
larger sampling. If the noise is locally Gaussian, this should be true independent of whatever 
technicality is involved in the multiplicative noise integration. Therefore, there would remain 
a residual current in the form of < F(x)/Y(x) > 7 ^ 0 as a result of having Boltzmann 
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distribution and some ad hoc cancellation of this current in this so called equilibrium will 
be required. 

In two previous papers [l2, 13] we have discussed these issues and it has been shown on 
the basis of the sole consideration of the non-existence of any average current over inhomoge¬ 
neous space in equilibrium that the stochastic force strength for equilibrium dynamics should 
be generalized to g(x) = r(x)^/2/csT/ < T(x) >. This choice also explicitly decouples the 
local damping T(:e) from the global diffusivity D = k B T / < T(x) > in the dynamics. 

With the above mentioned modified stochastic force strength which is a linear function 
of T(x) the over-damped dynamics (Langevin dynamics) becomes a stochastic problem with 
additive noise and the distribution one gets readily as 


P(x) = N exp 


< T(x) > 

k B T 


dx 


F(x) 

r(h 


(i) 


where N is normalization constant. It is important to note that, according to this present 
model, the local diffusivity does not characterize equilibrium. The local diffusivity can 
always be defined in relation with the local damping over a suitable local average, but, that 
does not necessarily mean that it has to determine the equilibrium of the system. It is the 
global diffusivity of the system defined over the whole finite inhomogeneous space of the 
system that features in the equilibrium dynamics. Now, using this distribution it’s straight 
forward to see that there is no average current in equilibrium and the distribution becomes 
the standard Boltzmann distribution when T(x) is a constant. There are previous attempts 
to generalize Boltzmann distribution based on the entropy proposed by Tsallis and in this 


regard the ref. 


16 


17] are interesting works involving nonlinear Fokker-Planck equation. 


In this paper we are going to show a full derivation of the generalized Langevin dynamics 
including the inertial term and the above mentioned modified noise strength starting from a 
Hamiltonian that includes bath degrees of freedom as well as the system. We will follow here 
a formalism due to Zwanzig Q . It will be shown that if there exists a local conformation 
dependent coupling of the system to the harmonic bath the mesoscopic dynamics of the 
system will have the stochastic noise of the modified form. In the following we first summa¬ 
rize the results of the generalized Langevin dynamics (which has been obtained in details 
in ref. 13|) in order to gain an idea about exactly what type of an equilibrium scenario is 
expected. After that, we would present a derivation of this generalized Langevin dynamics 
using the microscopic Hamiltonian. In the end we present a discussion. 
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I. THE GENERALIZED LANGEVIN DYNAMICS 


Consider the generalized Langevin dynamics including the inertial term as 
d 2 x 


m— = -mr(i)^ + F{x) + r(x)^ 


2mk B T 


■vit), 


( 2 ) 


< r(x) > 

where nr is the mass of the particle under consideration. In ref. [13|, the Fokker-Planck 
equation of the above model has been derived and the stationary distribution representing 
the equilibrium of the system is found out to be 

mv 2 < T(x) > 


P(x, v) = N\ 


1 m < r(x) > 1 

2nk B T r 3 / 2 (a:) 


exp 


< r(x) > 

k B T 


dx‘ 


r(z') 


exp 


2T{x)k B T 


(3) 


where v — -j and N is an overall normalization constant. The important thing to note here 
is that the velocity distribution is locally Gaussian with a position dependent width, which 
means that corresponding to this stationary distribution there exists a local temperature 
Tiocai = TF(x)/ < T(x) > which becomes the global temperature T on an average over 
the whole inhomogeneous space. For example, if the inhomogeneous space is the internal 
space of a protein molecule, its different conformations would see a conformation dependent 
temperature but, as a whole, the temperature of the molecule over macroscopic time scales 
is that of the bath which is T. The molecule would remain in equilibrium with the bath 
as it samples through its conformation space being driven by the equilibrium fluctuations 
which indicates the < r(x) > to be actually a weighted average which we would see to be 
indeed true as we look at equipartition of energy. 

Interesting to note that, although the conformations are characterized by local tem¬ 
peratures the average velocity at each conformation is zero because the velocity distribu¬ 
tion is everywhere explicitly locally Gaussian with a zero mean. In the derivation of the 


probability distribution in ref. 


13] the standard procedure of setting the damping and the 


stochastic part dependent probability current to zero has actually been followed. Going 
by Einstein’s definition of diffusivity, the local diffusivity in equilibrium of such a system 
is D(x) = ( fc <r(^)> )(f^)) = <fg> ’ which h identical to the global diffusivity. This is the 
most important signature of the present alternative theory over the existing ones. In the 
present context, equilibrium is characterised by this diffusivity and that entirely determines 
the stochastic nature of the problem. So, in equilibrium over an inhomogeneous space in¬ 
duced by coordinate dependent damping, it is not the temperature at all conformations 
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which is constant, rather it is the diffusivity which remains constant. This observation, if is 
experimentally proven, can potentially have a much bigger role to play in the interaction of 
complex molecules involving barrier overcoming processes in equilibrium. 


Following ref. 


13], let us have a look at the equipartition of energy and local power 


balance of the system. The local power balance in the presence of a conformation dependent 
temperature is indispensable in equilibrium. In the presence of the local power balance, 
there would remain no need for existence of any current to make unbalanced energy to flow 
between conformations (states) to maintain global energy conservation. In other words, the 
local power balance would ensure a local conservation of energy on noise average at each 
conformation in equilibrium. 

The equipartition can be easily seen to hold on average over the whole inhomogeneous 
space as 


< E ke > = 


mN m < T(x) > 


‘lixknT 


dx- 


dvv 2 exp 


mv 


T 3 / \x) 
2 < T(x) > 


exp 


2Y{x)k B T 


< Y(x) > 
k B T 
k B T 




T(x') 


x 


2 < T(x) > 


dxY(x)P(x) = —(4) 


where P(x) is the distribution one obtains by locally integrating over all velocities. Consis¬ 
tent with our previous assessment we indeed get to see that the system as a whole maintains 
a constant temperature T at which it can equilibrate with a heat bath. 

Note that, the computation of the equipartition implicitly presents the largest time scale 
(Macroscopic time scale) associated with the system over which the temperature seen by the 
system is the temperature of the heat bath. This is the time scale over which the system 
sees the whole of its allowed conformation space (inhomogeneous space) being sampled by 
its equilibrium fluctuations. Thermodynamically temperature is a quantity associated with 
a time scale below which one would see fluctuations in it. Equilibrium does not mean one 
would see the same temperature at all time scales. Compared to the above mentioned macro¬ 
scopic time scale, the mesoscopic time scale is the one over which the coordinate dependent 
damping, corresponding local temperature and the stochastic form of the noise are defined. 
This particular time scale is there for every stochastic model defined at mesoscopic scales. 
For a homogeneous stochastic systems this mesoscopic time scale captures the temperature 
of the bath because the space is same everywhere whereas it is the macroscopic time scale 
(depending upon the extent of the inhomogeneosu space) over which the temperature of the 
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bath is captured for an inhomogeneous stochastic system as considered here. 

The local power balance can also be seen to hold by computing the noise averaged input 
power locally as 

E in = l f rV ( X,t ^J ’ ( 5 ) 

where the v(t) is given by the dynamics and its stochastic part would only contribute to Ei n 
on a noise average. Thus, 


Ei, 


(2 k B TT 2 (x)\ 
(2 k B TT 2 {x)\ 

V < r(*) > ) 


e~ nx)t f dt'e Ex)t ' <v{tW) > 

Jo 

e" r(x)t [ dt'e Ex)t 'd(t - t ') 

Jo 


mT(x)v 2 (x) = E diss , 


( 6 ) 


where v 2 (x) is the mean square velocity at x. This balance is clearly overruling any energy 
current between conformations over noise average because the energy would remain locally 
conserved on noise average at each state although the states are characterized by local 
temperatures in equilibrium. The local temperature is an associated stochastic feature of 
the existence of the local damping when it exists as is taken in the model. 

Now, the question naturally arises is - does such a model have any microscopic basis? 
If it does have some microscopic basis then it should be derivable from some microscopic 
Hamiltonian where the above mentioned assumptions should reflect. The most important of 
all assumptions underlying the existence of such a model with coordinate dependent damping 
is the assumption that the system interacts with different amounts of bath degrees of freedom 
at different conformations (coordinates). In the following derivation of the model from 
microscopic Hamiltonian we will be using this pivotal assumption of differential coupling of 
the system with the bath. 


II. MICROSCOPIC DERIVATION OF THE MODEL 


Consider the Hamiltonian of the full microscopic system as 

V 2 . v- k 2 ( l(X)X \ : 


H = 


*(X) + ]T 


2 


Qi 


OJt 




+ 


Pt 


(7) 


We consider the coupling constant ^(X) is a function of the mesoscopic system’s coordinate 
X and 'j(X) is independent of the bath oscillator index i. The conformation of the system 
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is given by momentum V and coordinate X. 

Under this Hamiltonian the bath space is spanned by the bath degrees of freedom (har¬ 
monic oscillators) specified by coordinates g* and momenta p t . We want to find out the 
stochastic dynamics of the mesoscopic system which is harmonically coupled to the bath by 
integrating out g* and p*. The <f>(X) is a potential causing the presence of some conserva¬ 
tive force to the mesoscopic system, a 7 is the frequency of the ?'th bath oscillator where the 
masses of the bath and system are considered unity without any loss of generality. 

The equation of motion for the system turns out to be 

V 1 a$(X) 




( 8 ) 


/(X) /(X) ax 

where /(X) = — , Di(t ) = g,(t) — 7 t A '^j A 'h) anc [ an OV er-dot implies total time derivative. 

The equation of motion for the ith bath degree of freedom is 

Qi + ufa = q(X)X, (9) 

which can be solved to give 

, /N 7 \X(t)]X(t) . . , 7 [X(0)]X(0) Pi ( 0) 

Di(t) = qi(t )-—- = (g*( 0 )- j cosw ^ + -—~ smujit 


LOT 


UJ? 


UJi 


COS UJiit — s ) 


f{X)V{s)ds. 


( 10 ) 


Jo V 

Now, by replacing D^t) in the equation of motion of the system one gets the dynamics 
of the system as 

v i a$(x) rt 


/(X) /(X) ax 


k(t — s)V(s)ds + 


where 


and 


k(t - s) = f(X)J 2 


COSO Ji(t — s ) 


Wf 


</>(t) = V (qi(0) - 7[A(0) j X(0) ) cos ujjt + ^21 sin Uit . 

OJ; OJj 


( 11 ) 


( 12 ) 


(13) 


Let us first compute the quantity k(t — s) at X as 


k(t — s ) 


= lim / 

(-^max ^00 J q 

= lim 2 C 

ta max yoo 


doj 


g{u>,X)f{X) 


U! 


COSO j(t — s ) 


du cosa :{t — s ) 


0 


= lim 2C SiaUmax{ t s) =2 C6(t-s). 

<jJmax—yoo id — S 


t — s 

















In the above computation we have taken the density of the oscillators of frequency oj at X to 
be g(co,X) = which is consistent with the scenario that at different conformations the 
system interacts with a varied number of bath degrees of freedom to have a conformation 
dependent damping and temperature. This choice of g(oj,X) is a local generalization of 


the density of oscillators as should be taken to generalize Zwanzig’s 


homogeneous space (additive noise) to non-homogeneous space [l4j, [15]. For the additive 


ormulation for the 


noise (homogeneous space) case (i.e. 7 (A) = 7 a constant), the density of oscillators would 
be taken as g(cu) = 2C “ 2 in order to get a white noise following Zwanzig’s formulation. Note 
that when 7 (A) = 7 , we have f(X) = 7 which clearly demonstrates the correspondence. 

Following the approach of Zwanzig, let us consider Maxwell-Boltzmann (MB) distribution 
of the initial coordinates Di( 0) and momenta Pi(0) of the bath degrees of freedom. However, 
the bath degrees of freedom we are considering here are the ones interacting with the system 
which is at its initial (t = 0) position and momentum. When the system is in equilibrium 
with a state dependent temperature r[A(0)]T/ < r[A(0)] >, guessing that T(A") oc /(A) 
we take the temperature of the distribution of the bath degrees of freedom at t — 0 (as seen 
by the system in equilibrium with the bath) to be /[A(0)]T/ < /[A(0)] >. Here, essentially, 
what being considered is that the particular bath degrees of freedom interacting with the 
system at a given conformation over a mesoscopic time scale are the only ones seen by the 
mesoscopic system at that local temperature. So, an average over bath degrees of freedom 
at t — 0 should involve only these degrees of freedom which are seen by the mesoscopic 
system. Therefore, the statistics of the bath as seen by the system at the initial instant 
in the standard way [15] is specified by the correlations involving the local temperature 
is given by the correlations < p*(0) >=< _Dj(0) >= 0, < Pi{0)Pj{0) >= k <flx(o)}> 

< Di(0)Dj(0)u}iU}j >= fc </rx( 0 )] < A;(0)/Jj(0) >= 0. Using these averages at the initial 
instant it can be shown that 
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< c t>(t)(J)(s ) > = / ^D, ; ( 0 ) cosc Jit + ^—1 sincujt^ ^Dj(O) cosc OjS + ^ sin ujjf^j 

\ ^ij ^ 

V- k B f\X(0)]T s , 4 , . , . 

= — - Oij (cos Uy£ COS UJjS + Sill Uit Sill UJjS) 


hi 


■f < /[A"(0)] > UiUJj 


k B T 


y -COS LUi(t — S ) 


< fX( 0 ) > -T- u,; 

(k„T)(2C5(t - s)) k B Tkit — s) 


( 14 ) 


< f{X(0)\ > < !(X) > ' 

Using the expressions for the second moment of </>(£) and k(t — s ), one can finally write Eq.10 
in the form 


V(t) = - Cf(X)T(t) + Cf(X), 


2knT 




(15) 


dX ~ J v “ y ' v “ y ' v “ y \/ C < f(X) > 

where 77 (f) is a white noise of unit strength and the factor of two in the damping term is 
cancelled by the half appearing due to integration over half of the delta function S(t — s). 
Eq.15 is the equilibrium stochastic dynamics of a Brownian particle in an inhomogeneous 
space (the same as Eq.2) when we identify mr(Jf) = Cf(X). A multi-dimensional extension 
of the same should follow the same line of logic. 


III. EXISTING EXPERIMENT 

There is an experiment done by Faucheux and Libchaber jf5| to find out the variation 
of diffusivity of different Brownian particles moving under gravity between two confining 
planes. This experiment was based on direct visualization of trajectories of particles. In 
this experiment the whole 3-dimensional trajectory of a BP is projected on a horizontal 
x — y plane. The diffusivity thus obtained on the transverse plane is an averaged quantity 
along the vertical direction. This average transverse diffusivity on horizontal planes D B is 
then plotted against a parameter which is the average height of a particle of a given radius 
from the lower confining plate. Obviously, this parameter scale is equivalent to the scale of 
particle size and is an appropriate one to plot vertically averaged horizontal diffusivity. It 
was shown that the diffusivity of a particle varies by as much as 2/3 of its bulk value at 
the closest proximity of the lower confining plate under experimental conditions. Here, the 
proximity is in terms of the parameter mentioned above which is a reasonably good measure. 
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A nice thing about such a system is that the vertical (^-direction) viscosity rj(z) is ex¬ 
actly known and the transverse (x,y) components of the viscosities were also used upto the 
approximate order (r/ z) 5 where r is the radius of the BP. The derivation of the present the¬ 
oretical results completely rely on the conformation dependent coupling of the system with 
the bath degrees of freedom and no hydrodynamics of the bath giving rise to the variation 
in damping. However, the knowledge of the coordinate dependence of damping in this ex¬ 
perimental system makes it an attractive one to probe the present results. If hydrodynamic 
effects are not too strong, most probably, such a system can reveal experimental results to 
prove the present theory. 

This experimental result, shown as early as in 1994, in the opinion of the present author, 
has clear indication of limitation of the Boltzmann distribution in an inhomogeneous space. 
The experimentally obtained horizontal average diffusivity Dh as plotted against the average 
height of a particle should be a diffusivity obtained by averaging over the vertical ^-direction. 
Now, the analytic expression of the z dependence of the planar diffusivity D\\{z) being known, 
if Boltzmann distribution Pb(z) works, then Dh = Jq (1zD\\(z)Pb{z) where t is the vertical 
gap between two confining plates. It was clearly mentioned in this paper that this quantity 
does not match the experimentally observed average diffusivity. Then, by some particular 
double averaging (one local and the other global) the diffusivity is computed in this work 


which fits the experimental data (refer to the original paper [6] for details). 

In the opinion of the present author, if the D\\[z) is correct and the vertical distribution 
function is the Boltzmann distribution ( Pb(z )) there is no reason why the direct averaging 
should fail. The experimental diffusivity has been extracted by nicely averaging over many 
sections of equal temporal stretches of a trajectory projected on horizontal plane (x — y 
plane). This projection on x — y plane is where the direct average in the vertical direction 
gets counted. The whole trajectory considered in this experiment is large enough to be able 


to properly see the available space. The sma 
showing an almost straight line (Fig.3 of ref, 


1-time part of the mean square displacement 
indicates a very good averaging has been 


done. Now, this average cannot be different from a direct average of the planar diffusivity 
if D\\(z) is correct and the distribution is the Boltzmann distribution. This is a beautiful 
experiment, most probably, giving a clear indication of a breakdown of the Boltzmann 
distribution in an inhomogeneous space where the inhomogeneity is caused by the space 
dependent damping. 
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Let us try to understand why such an experiment can be done to directly probe the 
modified Boltzmann distribution as presented in this paper before we elaborate on a proposal 
of a similar experiment. In the present experiment, it is a single Brownian particle whose 
motion near a wall is tracked. There can be multiple reasons for the variation in the damping 
near a wall for such systems. One of the reasons could be the differential coupling to the 
bath degrees of freedom if there exists a gradient of bath degrees of freedom near the wall. 
Or, there can most probably exist other hydrodynamic reasons, sticking of the boundary 
layers of water to the lower glass plate etc. Whatever be the microscopic origin of this 
damping variation with average height near the wall, if the energy lost by the particle to 
damping is retained by the fluid (bath) without loss and that transfer is a reasonably quick 
and local process, then it is a good system to probe proposed results. By reasonably quick 
we mean here that the time scale of this transfer is much smaller than the mesoscopic time 
scale. Such a situation should follow the same modified statistics as shown in this paper. 


IV. PROPOSAL OF EXPERIMENT 

We propose here two experiments one of which can be done exactly in the same 
way as was done by Faucheux and Libchaber for a direct visualization of the trajec¬ 
tory. In such an experiment one can obtain the (stationary) probability distribution of 
the vertical coordinate z of the BP. Then one can fit this distribution to a modified 
Boltzmann distribution as proposed in this paper. The modified Boltzmann distribution 
Pmb{z) — N exp (—V(z)/D) is structurally the same as Boltzmann distribution with a dif¬ 
ferent potential V(z) = dz'([m — m w \g/Qnrri^z 1 )) where m is mass, r is radius of the 
particle under consideration, m w is the mass of water of same volume as that of the particle 
and 7](z) is the coordinate dependent viscosity for the motion in ^-direction. The normal¬ 
ization constant N should be found out for the whole span of the available space in the 
^-direction between confining plates. While fitting this distribution to the experimentally 
obtained one, the diffusivity D in the vertical direction can directly be measured in the 
experiment. The average global damping coefficient < T(z) >=< 6nrr](z) > will come out 
to be ksT/D where T is the actual (known) temperature of the bath. 

Note that, this diffusivity D = fcgT/ < r(z) > in the vertical direction is dependent 
on the size or mass of the particle or in other words is dependent on the average height of 
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the particle which has been considered as a parameter in the experiment of Faucheux and 
Libchaber. A direct measurement of the vertical diffusivity D of the system, for example, 
above and below the average height of the particle would be much more interesting result 
in such an experiment. If the present theory works here, one should get to see the vertical 
diffusivity remains the same above and below the average height of the particle although 
there is a variation in the damping. The important thing to keep in mind in the evaluation 
of this vertical diffusivity is that it is not the average diffusivity < KbT/Q urr/^z) > in this 
direction, rather, it is a diffusivity determined by dividing the kuT by < 67irr](z) >. It is 
the diffusivity corresponding to the average viscosity < r](z) > and is not the average of 
diffusivity in the vertical direction. 

The transverse planer diffusivity that Faucheux and Libchaber have plotted against the 
average height, according to the present theory, should be averaged using this new dis¬ 
tribution in the vertical direction. The viscosity on horizontal planes at a given height 
being uniform everywhere the conventional theory of Brownian motion would work on these 
planes. On the transverse plane the distribution is the Boltzmann distribution if any pla¬ 
nar confinement exists through a conservative potential. Otherwise, it is just conventional 
open diffusion on the transverse planes with diffusivity D\\(z). Now, the average horizontal 
diffusivity should be Du = J 0 * cIzD\\(z)PmB(z). If the direct average in this way matches 
the experimentally obtained one then that will also be the proof of existence of the modified 
distribution. 

A second experiment of different type where the present theory should be more applicable 
would be the one involving a complex molecule e.g. a protein or a star polymer immersed 
in a fluid at a constant temperature and staying far away from the boundary. A direct 
tracking of the conformations of the system over time should be done here. Then, the 
accessible surface area of the system to bath degrees of freedom can be estimated out from 
those conformations. This accessible surface area to bath degrees of freedom can constitute 
the coordinates of the one-dimensional space over which the diffusion of the system can be 
defined. This diffusion over the coordinate space of accessible surface area would then be 
the one-dimensional diffusion under the equilibrium fluctuations. 

One can directly probe here the constancy of diffusivity over the whole conformation 
space although the damping can change from conformation to conformation. These are 
systems probably where the damping variation with conformation is not as clearly known 
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as that in the experiment of Faucheux and Libchaber. However, if the diffusivity comes out 
to be a constant in such an experiment, then one can find out the effective potential from 
the fitting of the graph of probability distribution to the experimentally obtained one in this 
one-dimensional space. A comparison of the effective potential with the actual one can then 
give a conformation dependent viscosity back if that exists. The weighted average of this 
conformation dependent viscosity then should again give back the same constant diffusivity 
to complete the self-consistency check. 


V. DISCUSSION 

In the present paper we have considered a mesoscopic system in contact with a heat bath 
where the coupling of the system with the heat bath depends on the conformation of the 
system. Starting from this microscopic scenario of conformation dependent coupling of a 
mesoscopic system with a heat bath we have derived the mesoscopic dynamics of the system 
which is in the structure we have expected. The form of the mesoscopic dynamics actually 
differs from the other conventionally taken Langevin dynamics where the stochastic force 
strength is assumed as y/2mT(x)k,BT. 

We introduced the local MB distribution for the bath degrees of freedom just before 
Eq.14 and that is an assumption. It does not anyway mean that because of making that 
assumption the distribution of the mesoscopic system would also be a similar one. Rather the 
distribution of the mesoscopic system has to be found out once we know the dynamics of the 
system after integrating out the bath degrees of freedom. Since, the dynamics of the system 
comes out in the expected form corresponding to which we know the distribution we can 
write down the distribution of the whole system (comprising of the bath and the system) 
in a product form. Integrating out the bath part then (when its distribution is suitably 
normalized), we can always recover the distribution of the system from such a distribution 
of the whole system. Note the interesting limit where T(x) becomes independent of x the 
distribution for the system we actually get is standard MB for the system and the bath alike 
as is obvious over homogeneous space. 

The important point here to make is we cannot a priori assume the distribution to be 
standard MB for a system which has conformation dependent coupling to a bath. This has 
to be found out explicitly from the mesoscopic dynamics and the assumption of standard 
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MB distribution in equilibrium would be a wrong one to make. It works for the homogeneous 
damping does not mean that it should be a general procedure. Note that, considering the 
system and the bath together does not mean that one is at thermodynamic limit. The system 
goes to thermodynamic limit when its boundary goes to infinity keeping the densities finite. 
But, that is not the case here and we are within a mesoscopic domain. 

The generalized Langevin dynamics with the conventional noise strength mT(x)/csT 
gives a standard MB distribution for supposed equilibrium and it has all the merits like 
equipartition and local energy balance as shown here for our modified model.The reason 
for this correspondence most probably lies in the fact that, the conventional generalized 
Langevin dynamics can also be converted to an additive noise problem by rescaling the 
velocity by ^/T(x) which is done in ref. 13] by T(a:). This indicates that to fix the correct 
stochastic noise strength one has to take into account the over-damped dynamics limit. 
Here, the crucial characteristic feature that the conventional method lacks is that one cannot 
directly go to the over-damped dynamics by dropping the inertial term from the dynamics 
which, on the contrary, is straightforwardly doable with our model. 

In the conventional case one has to augment the dynamics that results from dropping the 
inertial term with other ad hoc terms to cancel current. There is absolutely no need for such 
things for our model. In fact, our model with the noise strength T(X) w ^ iere th e 
noise is a linear function of the damping is the only form which would result in an additive 
noise dynamics for equilibrium in the over-damped regime. This over-damped dynamics is 
absolutely free of any spurious current. All other noise forms including the standard one 
would remain multiplicative noise dynamics in the over-damped regime and that generality 
most probably indicates the common non-equilibrium character of them. 

We would like to conclude by mentioning that the knowledge of the exact equilibrium 
probability distribution of a class of systems is not only required for the systems in 
equilibrium, rather, is required for all the estimations in the linear response regime which 
includes a lot more phenomena. For example, if experimentally proved, the conformation 
dependent temperature of complex molecules under equilibrium fluctuations can influence 
their reactions which cannot be explained by the Boltzmann statistics. If a wrong distribu¬ 
tion is used for equilibrium the same gets involved in the weakly non-equilibrium regime as 
well giving rise to anomalies which should have not appeared with a proper theory. One 
of the long standing unsolved issue is Leviuthal’s paradox in protein folding. Possibly, a 


15 







modification of the equilibrium distribution of mesoscopic systems holds the clue as to how 
a path is cut over a rugged energy landscape to make the system find the global minimum 
corresponding to the native fold of a protein. 
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